Electric Dipole Polarizabilities of Hydrogen and Helium Isotopes 
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The electric dipole polarizabilities of '^H, ''He, and "^He are calculated directly using the 
Schrodinger equation with the latest generation of two- and three-nucleon interactions. These polar- 
izabilities are necessary in order to obtain accurate nuclear-polarization corrections for transitions 
involving S-waves in one- and two-electron atoms. Our results are compared to previous results, 
and it is shown that direct calculations of the electric polarizability of ^He using modern nuclear 
potentials are smaller than published values calculated using experimental photoabsorption data. 
The status of this topic is assessed in the context of precise measurements of transitions in one- and 
two-electron atoms. 

PACS numbers: 21.60.De; 21.45.-v; 21.10.Ky; 31.15.ac 



I. INTRODUCTION 



Theoretical calculations of transition frequencies in hy- 
drogenic atoms and ions have reached a level of precision 
where tiny corrections due to nuclear structure and dy- 
namics are necessary in order to interpret the results of 
high-precision measurements in these systems. This has 
largely been the result of recent improvements in Quan- 
tum Electrodynamics (QED) calculations [U, 0, S|- In 
many cases the experimental errors and estimated sizes 
of uncalculated QED corrections are much smaller than 
the nuclear corrections, and one can thus use those mea- 
surements (corrected for QED effects) as an experimen- 
tal determination of various nuclear quantities 0, [E| • We 
briefly review several such determinations. 

For S-wave hyperfine transitions in one-electron atoms 
and ions 0, @, Q , experimental precision is much greater 
than that of all theoretical calculations, while uncal- 
culated theoretical contributions to transition frequen- 
cies (including QED corrections) are significantly smaller 
than nuclear effects. The leading-order (i.e., largest) nu- 
clear contribution to these hyperfine transitions (called a 
Low moment Q) is determined by a correlation between 
the nuclear charge and current operators [1, 0|- Low 
moments may be further decomposed into Zemach mo- 
ments |9|] (viz., utilizing only ground-state expectation 
values of the charge and current operators) and polar- 
ization contributions (viz., including only virtual excited 
states between the two operators), both of which play 
significant roles. For the important proton (i.e., ^H) case 
the polarization effects are significantly smaller than the 
static (Zemach) corrections because the proton is much 
more difficult to excite than any nucleus [13, [HI, fl^ . 
Although exceptionally interesting, hyperfine transitions 
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are not the focus of this paper. 

The frequencies of transitions between S-states in hy- 
drogenic atoms and ions can be separated into a ref- 
erence value (essentially the Dirac transition frequency 
modified by reduced-mass effects) plus the much smaller 
Lamb shift. The Lamb shift contribution is dominated 
by QED corrections, but nuclear effects play a significant 
role in the best measured transitions. These nuclear cor- 
rections can be decomposed into nuclear finite-size cor- 
rections (i.e., determined by nuclear ground-state charge 
distributions) plus nuclear polarization corrections (viz., 
involving virtual excited states of the nucleus). The latter 
are typically dominated and determined by the electric 
polarizability, which reflects the distortion of the nuclear 
charge distribution as it is attracted by (and follows) the 
orbiting electron. 

The most accurate measurement of such a frequency 
was performed in Ref. [l3| for the 1S-2S transition in 
hydrogen, with a relative error of slightly more than 
1.4 parts in 10^'* and with an absolute error of 34 Hz. 
That error is slightly smaller than the estimated polar- 
ization correction of 60(11) Hz from Ref. JdJ, and is 
much smaller than the size correction of about 1000 kHz. 
The mismatch in the sizes of these nuclear corrections 
reflects both their different dimensional structure [l| and 
the fact that the proton is difficult to excite (compared 
to a nucleus), even though the proton size is not signif- 
icantly smaller than that of light nuclei. If one turns 
the problem around and extracts the proton-size correc- 
tion from the experimental transition frequency [l|, one 
obtains a value for the proton r.m.s. charge radius of 
{r'^yjij^ = 0.877(7) fm, which agrees with a recent di- 
rect determination of that quantity from elastic electron- 
scattering data [H, [3: {r^)li^ = 0.897(18) fm. Both 
the polarization-correction and experimental errors are 
much smaller than the Rydberg constant error [l| , which 
dominates the uncertainty in the hydrogen atom analysis. 

A similar analysis of transitions from the 2S state 
in deuterium to a variety of S and D states [l| leads 
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to a value for the deuteron charge radius of (r^)J{^ = 
2.1402(28) fm, which is consistent with the electron scat- 
tering value [13, m of (r2)^/^ = 2.130(10) fm. We 
note that this is the full charge radius (containing the 
finite sizes of the proton and neutron constituents), and 
that the atomic value has an uncertainty nearly 4 times 
smaller than the value obtained directly from electron 
scattering. 

The determination of the difference in transition fre- 
quencies between hydrogen and deuterium for identical 
transitions can be used to test our understanding of small 
contributions to the charge radius of the deuteron 
In such a difference nuclear-mass-independent terms (in- 
cluding difhcult-to-calculate QED contributions) cancel, 
which greatly simplifies the analysis. Because the finite 
size of the proton contributes linearly to the deuteron 
mean-square radius (which is the nuclear quantity that 
determines the dominant nuclear-size correction in an 
atom), it largely cancels out in the frequency differ- 
ence. Higher-order proton-size corrections and neutron- 
size corrections are relatively small and tractable. The 
transition-frequency difference (dominated by calculable 
reduced-mass effects) was reported in Ref. [13| for 1S-2S 
transitions with a relative error of 2.2 parts in 10^° and 
an absolute error of 0.15 kHz. The nuclear electric po- 
larizability of deuterium contributes 19.26(6) kHz [20], 
which is more than two orders of magnitude greater than 
the experimental error, while the deuteron-size correc- 
tion is greater than 5000 kHz. The weak binding of the 
deuteron makes possible the calculation of the bulk of 
the polarization and nuclear-size corrections in terms of 
a few well-measured parameters. The tiny remaining size 
correction includes statistically significant contributions 
to the nuclear charge radius arising from meson-exchange 
currents and relativistic corrections [3, [HI , which are un- 
obtainable from other types of experiments such as elec- 
tron scattering. Obtaining this sensitivity to fine details 
of nuclear dynamics depends on accurate estimates of the 
deuteron electric polarizability. 

Measurements of S-wave transition frequencies in '^H, 
•^He, and ^He atoms are not yet as accurate as those 
described above, nor are the necessary theoretical calcu- 
lations for He atoms. It may be possible to improve [l^] 
both to the point where nuclear physics information can 
be extracted, particularly information about the r.m.s. 
charge radii. As reviewed and updated in Ref. [2^ . 
on the other hand, isotopic differences in transition fre- 
quencies for helium and singly ionized lithium isotopes 
now have the required experimental and theoretical sen- 
sitivity. The latter sensitivity is greatly enhanced by 
the cancellation of nuclear-mass-independent relativistic 
and QED corrections in isotopic differences. In com- 
plete analogy to the hydrogen-deuterium case, calculable 
reduced-mass effects dominate the frequency differences, 
leaving nuclear contributions as the residue after rela- 
tivistic and QED contributions are subtracted. There has 
been considerable recent interest in the isotope shifts of 
^He [2l[2|, ^He [2i,[23,[23, and «He [H transitions rel- 



ative to those of ^He. In each case the value of the r.m.s. 
charge radius of that isotope has been extracted relative 
to the charge radius of ''He [2^,[33l- The nuclear polariz- 
ability correction to the "^He - ^He isotope-shift frequency 
(the best measured of the He isotope shifts ) is about 2/3 
of the 3 kHz experimental uncertainty [25, [2^ . |2^ , while 
presently only a marginal influence |23| on the others, 
but future improvements should require reliable values 
of their electric polarizabilities (as was the case for the 
deuteron), and that is the purpose of this paper. 



II. CALCULATIONAL TECHNIQUES 

The electric (dipole) polarizability of a nucleus (or 
atom), ckE, is defined by 



\{N\DM\' 
En — Eq 



(1) 



where a is the fine-structure constant, Eq is the energy of 
the ground-state |0), En is the energy of the Nth excited 
state, \N) (all of which are in the continuum for few- 
nucleon systems), and is the component of the (non- 
relativistic, in our case) electric-dipole operator in the z 
direction, which generates the transition between those 
states. The definition (1) can be rearranged into the form 
of a sum rule 
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where cr"'^(a;) is the nuclear cross section for photoabsorp- 
tion of unretarded-dipole (long- wavelength) photons with 
energy uj, and wth is the threshold energy for photoab- 
sorption. The inverse-energy weightings in Eqns. ([1]) and 
([2]) lead to significant sensitivity of as to the threshold 
energy, wth, which depends on nuclear binding energies. 

In order to obtain the nuclear energy spectra and the 
wave functions involved in the calculation of the elec- 
tric polarizability (Eqn. ([1])), we use the no-core shell 
model (NCSM) in relative coordinates [3l| to solve the 
Schrodinger equation. The NCSM is a flexible approach 
to solving the few- and many-nucleon problems, and it 
has been extensively used in studies of s- and p-shell nu- 
clei [si, [H H [IS Hi]. In the NCSM the nuclear wave 
functions are obtained by the diagonalization of an effec- 
tive Hamiltonian in a finite basis constructed from har- 
monic oscillator (HO) wave functions. The truncation 
of the model space is taken into account via an effec- 
tive interaction derived by means of a unitary transfor- 
mation. Either local or non-local nucleon-nucleon (NN) 
and three-nucleon (NNN) interactions can be used in 
the Schrodinger equation. The effective interaction is 
constructed in a cluster approximation, which must be 
truncated for practical reasons. The truncation of the 
model space is determined and labeled by the number 
of excitations, Nmax, above the non-interacting state. 
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We test convergence by plotting calculated quantities vs. 
Nmax, and those quantities should approach their cor- 
rect asymptotic values as Nmax becomes infinite. Thus 
by observing the convergence of observables as a function 
of Nmax, we can determine their values. 

In this paper we compute the "^H, "^He, and '*He electric 
dipole polarizabilities starting from a nuclear Hamilto- 
nian derived within the framework of (QCD-based) Chi- 
ral Perturbation Theory (including the Coulomb interac- 
tion between the protons) . We adopt the nucleon-nucleon 
interaction at next-to-next-to-next-to-leading order (or 
N'^LO) of Ref. and the three-nucleon interaction at 
next-to-next-to-leading order (or N^LO) [sl, [1^ in the 
local form of Ref. [3.f!|. The accuracy of these nuclear 
interactions for s- and p-shell nuclei was investigated ex- 
tensively in the same NCSM framework in Ref. [13] ■ In 
particular, the experimental binding energies of '^H and 
•^He are reproduced with high accuracy (viz., within 8 
keV, or about one part per thousand) [H, [H, , while 
the '*He ground-state energy is within a few hundred 
keV of experiment (i.e., within approximately 1%). This 
residual discrepancy with experiment reflects our current 
uncertainty on the underlying nuclear dynamics. These 
modern nuclear forces therefore provide an accurate de- 
scription of the structure of the nuclides considered here 
(•^H, ■^He, and ^He) as well as the total photoabsorption 
cross section of "'He [s^ (discussed in the next section). 

For each nucleus we first solve the few-nucleon 
Schrodinger equation in order to obtain the ground-state 
wave function, which can be accurately described in a 
large HO basis. We next rearrange Eqn. ^ according to 
Podolsky's technique [4ll| , which allows the ground state 
to be used as the driving term for the Lanczos-monient 
method [13, , which is our numerical method of choice 
for solving the Schrodinger equation. This trick allows 
us to work only with bound-state quantities and to by- 
pass the much more difficult approach of computing a 
response in the continuum. A detailed description of this 
method in a NCSM framework was given in Ref. [i^l . 



III. RESULTS OF CALCULATIONS 

In Figs. [T]-[31 we show the running of the electric po- 
larizability with the truncation parameter for the model 
space, Nmax- Different HO frequency parameters, 11, re- 
sult in different convergence patterns for the electric po- 
larizability, and this fact is especially visible for small 
^max values. As shown in Figs. [IHSl results obtained 
using smaller HO frequencies (equivalent to larger-length 
oscillator parameters, defined hy b — l/vn7/vT2) converge 
faster. Long-range operators (such as the electric dipole 
operator) are thus better described using smaller values 
of O in the smaller model spaces. Moreover, better sam- 
pling of the low-lying excited states (the most important 
states for the calculation of the electric polarizability) 
is obtained for small values of VI. Although not shown, 
other operators converge faster for larger HO frequencies. 



Binding energies, for example, achieve the fastest conver- 
gence for a HO length parameter h on the order of the size 
of the nucleus considered. For each observable the results 
that arc obtained with different values of nevertheless 
approach a single asymptotic value for large Nmax- Un- 
certainties in determining that asymptotic value lead to 
error estimates in Table HI 

The upper panels of Figs. [T]-[3] present results obtained 
by neglecting three-nucleon interactions, while results 
that include three-nucleon interactions are shown in the 
lower panels. We note, however, that because binding 
energies are significantly smaller than experiment in the 
absence of three-nucleon interactions, the values obtained 
with only NN interactions are about 10-25% larger than 
the results obtained when NNN interactions are included. 
This is partly the effect of having too small a value for 
wth, which emphasizes smaller values of the energy de- 
nominators in Eqns.lH]) and 

The stronger binding of '^He than that of the three- 
nucleon systems and our slightly less accurate reproduc- 
tion of the ^He experimental binding energy may affect 
our results for aE. In order to probe this possibility we 
have performed three calculations with slightly different 
parameter sets in the three-nucleon force. The specific 
results are given in [ist . 

We expect from Eqn. ^ that cke should scale di- 
mensionally like the square of the nuclear size divided 
by the binding energy. Moreover the square of the size 
should scale roughly like the inverse of the binding en- 
ergy (this statement is highly accurate for the weakly 
bound deuteron) . The resulting dependence on the bind- 
ing energy should be roughly like the inverse square, and 
our slight over-binding could diminish Q!e by as much 
as 2%. Our calculation that neglects the three-nucleon 
force in ^He results in a value of q;e — 0.0822(5) fm'^ 
(22% higher than a result of 0.0673(5) fm'^ that incorpo- 
rates this force) , while reducing the binding energy from 
28.50(3) to 25.39(1) MeV (an 11% decrease). A similar 
effect is also seen in the '^He and ^^H calculations. 

A more difhcult problem is that stronger binding 
emphasizes nuclear corrections of relativistic order, in- 
cluding corrections to the electric-dipole operator from 
meson-exchange currents, which are determined by de- 
tails of how the nuclear forces are constructed [46;] . This 
effect could be as large as ±2% for the well-bound *He, 
but has been little studied and takes us far beyond the 
scope of this work. We will incorporate these uncertainties 
into our results in the next section. 



IV. COMPARISON WITH OTHER WORK 

Our calculations of the electric polarizabilities of three- 
and four-nucleon isotopes of hydrogen and helium are 
summarized in Table [H together with those of others us- 
ing different two-nucleon and three-nucleon forces. We 
have restricted our own entries to those that incorpo- 
rate three-nucleon forces and hence have accurate bind- 
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FIG. 1: [Color online] The dependence of (the electric po- 
larizability) oe of ''H (in units of fm'') on the model-space 
truncation parameter, Nmax ■ The results have been obtained 
using (a) NN interactions only, and (b) NN+NNN interac- 
tions. Each curve is obtained using a different frequency pa- 
rameter for the basis states, shown in the legend in MeV. For 
sufficiently large Nmax each result should be independent of 
that frequency. 



ing energies, especially for the three-nucleon systems and 
slightly less so for "^He. For completeness in the table we 
have also included the deuterium and ^He cases, which 
were not treated in this work. 

Only one other calculation of for '^H exists [l^l, and 
our result is in agreement with that calculation. 

Calculations for the electric polarizability of '^He (isl . 
[s^ l are in agreement within their uncertainties, and are in 
reasonable agreement with the determination of Ref . [5^ , 
but not with Ref. [52] . We note that if charge symmetry 
were exact in the three-nucleon systems, the Hamiltoni- 
ans and polarizabilities of '^H and '^He would be identical. 
Under the charge-symmetry operation that relates the 
two nuclei the dipole operators in Eqn. ([T]) would each 
develop a minus sign, while the radial wave functions 
and Green's functions would be identical. Most of the 
charge-symmetry violation in these systems is caused by 
the repulsive Coulomb interaction between the two pro- 
tons in ^He. We note that our uncertainties for these 
two nuclei are also different. The repulsive Coulomb in- 
teraction in ■^He leads to a larger radius for that nucleus, 
and that may be the source of the larger uncertainty. 
Matrix elements of infrared operators (i.e., those like the 
dipole operator that are most sensitive to large distances 
from the center of a nucleus) converge more slowly in 
the NCSM than do short-range operators, which can be 
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FIG. 2: [Color online] Same as in Fig. \T\ but for ^He. 
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FIG. 3: [Color online] Same as in Fig. [T] but for ^He. 

successfully renormalized [13, HI] . 

The uncertainty in the underlying nuclear dynamics 
(rather than the uncertainties reflected in the conver- 
gence plots) dominates the error estimate of our cal- 
culated electric polarizability of ^He (see (1^). After 
correcting our three values (corresponding to three nu- 
clear force models) for overbinding [isj we average them 
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TABLE I: Values of the electric polarizability of light nuclei, 
both theoretical and experimental, in units of fm'^. The exper- 
imental results have been determined by nuclear experiments, 
including the use of experimental photoabsorption data in 
Eqn. ([2]). No uncertainties were given for the ^H, '^He, and 
*He calculations in [47l. |48||. but they are likely to be smaller 
than about 10%. The ®He result is a hybrid calculation re- 
lying on some theoretical input and we add it here for com- 
pleteness. Results from the present calculation have no listed 
reference. Our three separate results ^45y] for *He have been 
combined in the table and are discussed in the text near the 
end of this section. The result of Ref. [4^ for the deuteron is 
an Effective Field Theory calculation. The errors for the two 
deuteron calculations are not independent and should not be 
combined. 



Nucleus 


aif'" (fm^^) 


ref. 


cxp /r 3\ 


ref. 


2h 


0.6328(17) 
0.6314(19) 


[20] 
[49] 


0.61(4) 
0.70(5) 


[50] 
[51] 




0.139(2) 
0.139 


[47] 






^He 


0.149(5) 

0.145 

0.153(15) 


[48] 

M 


0.250(40) 
0.130(13) 


[52] 
[53] 


*He 


0.0683(8)(14) 

0.0655(4) 

0.076 


[56] 
[48] 


0.072(4) 
0.076(8) 


[55] 

Mi 


^He 






1.99(40) 


[54] 


and use 


their spread 


as our 


direct uncertainty, 


with 



an additional 2% systematic uncertainty from the nu- 
clear dynamics (discussed in Section III). This produces 
aE = 0.0683(8)(14) fm^, which is listed in Table D Our 
result there is significantly smaller than most of the cor- 
responding results, although just at the limit of the es- 
timated uncertainties. We are, however, in fairly good 
agreement with a recent calculation by Gazit et al. [5a |. 
which predicts a slightly smaller polarizability, but also 
corresponds to a slightly overbound model. We note that 
Ref. [48| used a very primitive nuclear force model and 
that those results are superseded by those of Ref. [56j . 
References [s^l and [lH used fits to experimental pho- 
toabsorption data and Eqn. ([2]) in order to obtain their 
results. Values obtained from a direct solution of the 
Schrodinger equation are therefore at some variance with 
those calculated using experimental "^He photoabsorption 
data. 

Measurements of ''He photoabsorption in the near- 
threshold region have been controversial over the years, 
particularly with respect to the height of the cross sec- 
tion at the peak, for which one can find differences of up 
to a factor of two between different experiments (e.g., see 



Ref. [36| and references therein). This makes it very chal- 
lenging to extract an accurate and unambiguous value 
of the ^He electric polarizability from the measured "^He 
photoabsorption cross section using Eqn. In contrast 
there has been substantial recent progress in theoreti- 
cal calculations of the ''^He photoabsorption cross section. 
Predictions obtained using high precision NN and NNN 
interaction models (including the ones used in this work) 
all lie in a rather constrained band [36j . in remarkable 
contrast to the large discrepancies present among differ- 
ent experimental data. This gives us confidence that our 
prediction for the "^He electric polarizability, obtained by 
direct solution of the Schrodinger equation, will prove 
to be more accurate than those obtained using existing 
photoabsorption data. 



V. CONCLUSION 

We have used the latest generation NN and NNN in- 
teractions in a NCSM framework in order to obtain accu- 
rate three- and four-nucleon solutions of the Schrodinger 
equation. Using the Lanczos-moment method to imple- 
ment Podolsky's technique [JTl for treating second-order 
perturbation theory, we have calculated the electric po- 
larizabilities of "^H, "^He, and '*He. Our result for "^H is 
in excellent agreement with that of Ref. [i^l , while that 
for '^He is in good agreement with previous work. Our 
best estimate of 0.0683(8) (14) fm^ for -^He based on di- 
rect solutions of the Schrodinger equation is at the lowest 
end of the calculations that used experimental photoab- 
sorption data directly in Eqn. ([2]) . Future calculations for 
other light nuclei such as ^He and ^Li should be tractable, 
but would require a change of basis for the NCSM. For 
nuclei with mass numbers greater than five, a Slater De- 
terminant basis is much more efficient than the relative 
coordinate approach used in this work. 
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